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ABSTRACT 

Fokker-Planck models are used to give estimates for the retention fractions for newly- 
born neutron stars in globular clusters as a function of kick velocity. These can be 
used to calculate the present day numbers of neutron stars in globular clusters and 
in addressing questions such as the origin of millisecond pu lsars. As an example, the 
Population I kick velocity distribution of Lyne & Lorimer ( 1994 ) is used to estimate 
the retained fractions of neutron stars originating as single stars and in binary systems. 
For plausible initial conditions fewer than 4% of single neutron stars are retained. The 
retention fractions from binary systems can be 2 to 5 times higher. The dominant 
source of retained neutron stars is found to be through binary systems which remain 
bound after the first supernova, ie. high-mass X-ray binaries. The retained fraction de- 
creases with an increasing number of progenitors, but the retention fraction decreases 
more slowly than the number of progenitors increases. On balance, more progenitors 
give more neutron stars in the cluster. 

Key words: globular clusters: general - stars: neutron 
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1 INTRODUCTION 

Globular clusters are now known to contain a population 
of neutron stars. There are 34 pulsars observe d in globula r 
clusters, 30 of them of the millisecond variet y (|Lyne 1995 ), 
and 12 low-mass x-ray binaries (LMXBs) ( prindlay 199^ ) 
which probably contain neutron stars. 

The generally accepted scenario for neutron star forma- 
tion and evolution gives their origin in Type II supernovae 
and for them to slowly spin down. The millisecond pulsars 
are thought to be the result of recycling by accretion of mat- 
ter from binary companions. For globular clusters, as well 
is in the Galaxy at large, this picture has been the matter 
of some controversy. The presumed progenitors of the mil- 
lisecond pulsars are the LMXBs, but some (see for example 
Bailyn & Grindlay 1990) feel that this is inconsistent with 
the relative numbers of the two objects. Rather than being 
spun-up neutron stars, the millisecond pulsars would, in this 
picture, come from the accretion induced collapse of a white 
dwarf pushed over the Chandrasekhar limit. (See the review 
of Bhattacharya & van den Heuvel 1991 for a thorough dis- 
cussion of these issues.) One of the key issues in this debate 
is the rate at which stellar interactions cause neutron stars 
to end up in binary systems leading to mass transfer and 
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being spun up. An essential ingredient in estimating this is 
the total number of neutron stars present in the system. 



Lyne & Lorimer (1994) have recently reevaluated the 
distribution of space velocities of young field pulsars. Based 
on their sample the mean birth, or "kick" , velocity of neu- 
tron stars is 450 km s" 1 . The typical escape velocity of a 
globular cluster is of order 20 km s~ , so if the globular clus- 
ter neutron stars are born with similar velocities, we would 
expect very few to remain in the clusters. 

Verbunt & Hut (1987) attempted to estimate a global 
retention fraction for globular cluster neutron stars. They 
assumed a typical escape velocity of 30 km s _1 and used the 
l ower kick velocity distribution of Lyne, Anderson, & Salter 
HH). They estimated the retention fraction be ~ 0.15. 
This t reatment was updated by Hut, Murphy, & Verbunt 
(1991) who used a collision- weighted average escape velocity 
of 47.7 km s _1 and estimated the retention fraction to be 0.3. 
These estimates are unsatisfactory since they are based on 
a single escape velocity and do not account properly for the 
variation in escape velocity within each cluster and between 
clusters. There is a need for a more realistic treatment of 
neutron star retention using dynamical models of globular 
clusters. This is all the more necessary in view of the increase 
in the mean kick velocity in Lyne & Lorimer (1994). 

The aim of this paper is to provide curves giving re- 
tention fractions as a function of kick velocity for a series 
of dynamical Fokker-Planck models. These are presented in 
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the next section. As an example of their use, in §^ I will cal- 
culate the overall re tentio n rates for single stars assuming 
the Lyne & Lorimer (1994) velocity distribution, and the re- 
tention rates of bi nary f ragments using the results of Brandt 
& Podsiadlowski fll995|). 



2 RETENTION FRACTIONS 

Consider a star with velocity v at a position in the cluster 
where the escape velocity is v e which receives a velocity 
kick Vk- We assume that both velocities are isotropically 
distributed so that the angle 9 between them is distributed 



P{0) = ~sin0. 

The resultant velocity, v r , is 

v r — \J v 2 + w? + 2vvk cos 8. 



(1) 



(2) 



The star will remain bound to the cluster provided that 
v r < v e , ie. so long as 



COS# < COS0 e 



2wk 



(3) 



The probability that the star will be retained by the cluster 



P(iet\v,Vk,v e ) = 



P{G)d6 



0. 



— 



1 < 

-i| < — ; 
l > -2-. 



(4) 



Note that Vk and ?; in eq. (^) can be written in terms of v e , 
so let fi = w/t; e and -fife = Vk/v e and then 



P(ret|fi,fi fc ) 



1, 

l-(v~v k ) 2 
4vv k 

0, 



fi fc — 1 < —v; 

\vk — 1| < fi; 
fife — 1 > fi- 



(5) 



If fife > 1, then all stars with v < Vk — v e will escape from 
the cluster, while if fife < 1, then all stars with v < v e — Vk 
are retained. All stars escape if fife > 2. 

In order to get the total retention probability as a func- 
tion of Vk for a given star system, we need the number of 
stars as a function both of position and velocity. For the pur- 
poses of this calculation, I assume a spherically symmetric 
cluster with an isotropic velocity dispersion. The distribu- 
tion function is then a function of the energy alone. The 
density and potential can be derived from the distribution 
function, as can the velocity distribution at each radius. 

In general, for a known distribution function, f(r,v), 
the number of stars retained at radius rasa function of fife 
is 



(6) 



N Tet {r,Vk) = 47T / P(ret.\v,v k )v f(r, v)dxi. 
Jo 

Using eq.(||), 

N Ict (r,v k ) = H(l-v k )iTV [ "v 2 f(r,v)dv 

Jo 



Vk 



v(l-(v-v k f)f(r,v)dv, (7) 



where H(x) is the Heaviside unit function. For a distribution 
function in energy f(E), as is appropriate for the isotropic, 
spherically symmetric models to be discussed below, we can 
rewrite the number retained as 



N le t(r,v k ) = 



47T 



(2 |#r) | 



ff(l-fih) / (1 - 2A)2 f (A) dA 



+- 



(l-2A)5/(A)dA 



The energy has been normalized by the local escape velocity 
so E = 2<f>(r)A, where the potential goes to zero at infinity. 
The limit on the integration is given by I = fife (2 — Vk)/2 and 
is a function of radius since fife depends on the local escape 
velocity. Equation ^ can then be integrated over all radii 
and divided by the total number of stars to give the total 
fraction of stars retained as a function of v k , /rot(ffc). 

The approach taken here was to evaluate the total frac- 
tion of stars retained using Monte Carlo techniques. To be- 
gin with , I consider the isotropic, single-mass, Michie-King 
models ( King 19*6^ ) . These form a single parameter family 
in terms of their degree of central concentration. The pa- 
rameter is often given as the dimensionless central potential 
Wo- For a series of models in this family, the density distri- 
bution and potential were calculated. For each test value of 
fife (the models are dimensionless and the local escape veloc- 
ity scales with the central escape velocity v ec ), stellar posi- 
tions were drawn at random with a distribution consistent 
with the radial number density. For each position, the veloc- 
ity distribution was computed from the lowered-Maxwellian 
distribution function. The relative angle between the veloc- 
ity vectors was chosen at random from the distribution in 
eq. ([!]) and the v r was computed. The retention fraction for 
the chosen fife is then the ratio of test stars with v r < v e (r) 
to the total number of test stars. Sufficiently large samples 
were used to ensure numerical precision. 

Figure [l] shows the retention fractions as a function 
of Vk /ve c , v £c is the central escape velocity, for a range of 
Michie-King models. The retention fractions are tabulated 
in Table |l] Note that the lower concentration models are 
most able to retain their neutron stars for a given Vk/v e< . 
since they have a greater number of stars within their cores. 
On the other hand, a low concentration model must have a 
much higher mass than a high concentration mo del t o have 
the same central escape velocity. This scales as yj Mjr where 
M is the total mass of the system and r is some scale radius. 

In order to make more realistic estimates of the reten- 
tion fractions, I have done similar Monte-Carlo integrations 
of eq. (0) for a series of Fokker-Planck models. These mod- 



els are as described in Drukier (1995) with the differences 
discussed here. The models are based on direct numerical 
integration of the orbit-averaged form of the Fokk er-Planck 
equation following the technique of Cohn ( 1980 ). A mass 
spectrum is included as are the effects of stellar evolution. 



(8) 
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Since only the early time evolution of the models is re- 
quired, the tidal stripping and binary heating as described 
in Drukier ( |l995j ) have been turned off. Stellar evolution 
has been handled by having each bin evolve over a time 
interval corresponding to the lifetimes of stars at the bin 
boundaries. Since even the most massive stars take time to 
complete their evolution, the earliest stages of the evolution 
of the models are governed by two-body relaxation alone. 
Dynamical friction can lead to a large degree of initial mass 
segregation when stellar mass-loss and neutron star forma- 
tion commences. 

The treatment of stellar evolution has been modified 



from that described in Drukier (1995) by adding an auxil- 
iary mass species to hold the distribution function for the 
remnants as they evolve. At each time step an appropriate 
fraction of the progenitor distribution function is transferred 
to the remnants. The mean mass of the two bins still changes 
linearly over the time interval for evolution, but each of the 
two populations of stars, with their very different masses, are 
free to evolve separately. At the same time as the transfer is 
done at each time step, the transferred distribution function 
and its associated density profile are saved, together with 
the current potential. These are then used in the integra- 
tion of eq. (ph at that time step. An appropriately weighted 
sum is computed to give the overall retention fraction for 
the model cluster. 

The mass range adopted for these models is from 32 
M0to 0.1 Mq. Neutron star progenitors are assumed to have 
masses above 8 Mq. It is likely that in Population II stars 
with masses as low as 5 MqwUI give rise to supernovae, but 



whethe r th ey leave a remnant or not is unknown (Weide- 



mann 1993). For such a large mass ratio, dynamical friction 



can be an important effect. I have used four logarithmically 
spaced bins for the neutron stars in these models. Tests us- 
ing only one or two bins gave substantially the same results, 
as did models with a lower mass limit of O.15M0. 

The models are assumed to start as Michie-King models 
with Wo = 3, 5, or 7; initial mass M = 10 5 , 5 x 10 5 , 10 6 , or 
5 x 10 6 Mq; and limiting radii r; between 15 and 65 pc with 
one model at 97 pc. The grid was not fully sampled. The 
initial mass function (IMF) is assumed to take the form of 
a power-law 



N(m) dm tx m 



dm, 



(9) 



with mass-spectral-index x; x = 1.35 for the Salpeter mass 
function. Models with x = 1 and x — 2 were run. For x > 2 
there are very few massive stars in the cluster to begin with. 

The initial parameters of the models run are listed in 
Table |^ and Table ^ sorted by initial half-mass relaxation 
time. The tables give the run identification number; the 
mode l parameters Wo, ri,x , and M; the half-mass relaxation 
time (Spitzer & Hart 1981); and the number of neutron-star 



progenitors. Only the models with relaxation times less than 
20 Gyr have been included. For the models listed in Table ^| 
the retention fractions as a function of Vk are listed in Ta- 
ble ^| The table columns are arranged by model number, 
and give the retention fractions for a range of velocities. The 
corresponding velocity in km s~ is the index number of the 
retention fraction multiplied by the number under the model 
number. I have listed the initial parameters of the models in 
Table ^| for completeness. These models retain none of their 
neutron stars. Their escape velocities are less than about 8 



km s _1 so that all neutron stars born with velocities greater 
than about 15 km s _1 escape. 

The effects on the retention fraction of varying the 
model parameters are shown in Fig. ^| At a given value 
of Vk the fraction retained increases with decreasing cluster 
size, increasing mass and increasing concentration. Models 
with smaller sizes, larger masses, or higher concentrations 
have larger escape velocities and can thus retain a higher 
fraction of stars at a given velocity. The variation with r; 
and M are comparable, but M has a much larger range in 
globular clusters. The variation with concentration is much 
the same as for the Michie-King models in Figure ^ when 
proper attention is paid to the scaling. 

For all values of Vk, the models with x — 2 have a higher 
f re t than the corresponding models with x = 1. Since the 
x — 2 models have fewer massive stars than the x = 1 
models, they lose less mass during the supernovae. Hence, 
the mean escape velocity over the period of massive star 
evolution is higher for the x — 2 models and more neutron 
stars are retained. The difference with x is fairly small and 
the additional retention fraction is tightly connected with a 
lower number of neutron-star progenitors. The net result is 
that the total number of neutron stars in an x — 2 cluster 
will be lower than in an x — 1 cluster. 

Although it is not obvious from Figure ^, apart from a 
velocity scale factor the retention curves for each x and Wo 
are very similar. Further, they are quite similar to curves for 
the appropriate Wo in Figure |l| To test this, I have com- 
pared the retention fraction curves with the curves for the 
appropriate Michie-King models. I have fit a scale velocity 
between the two curves and then computed the absolute 
maximum difference in f ret . For x = 1 the maximum differ- 
ence is 0.03 for model 1 and for x = 2 the maximum differ- 
ence is 0.05 for model 3. The size of the maximum difference 
generally decreases with increasing t T h, but the relationship 
is a complicated combination of the structural parameters 
Wo, n, and M. 

On the other hand, we can estimate the scaling velocity 
from the parameters as 



Ve(W )y/f M M/n, 



(10) 



where /m is the fraction of the total mass remaining after 
the end of neutron star formation and is 0.776 for x = 1 and 
0.991 for a; = 2. Values for v e (W o) are given in Table || for M 
in Mq and n in pc. (The King (1966) concentration param- 
eter c = log(r;/r c ) has also been listed for those unfamiliar 
with Wo-) This value of the scale velocity is within a few per- 
cent of the best fit estimate, but decreases with respect to 
the best fit value as the relaxation time decreases. Another 
way to look at this is that for short relaxation times, re- 
tention fractions calculated using the parameter estimate of 
the scale velocity and the Michie-King retention curves are 
systematically smaller than the retention fractions from the 
Fokker-Planck models themselves. This can be attributed to 
dynamical friction which enhances the number of neutron- 
star progenitors in regions with large escape velocities. The 
limited amount of two-body relaxation which takes place 
before stellar evolution mass loss begins also enhances the 
central potential and the neutron star retention rate. There 
is a small difference in the shape of the curves for x — 2 
and Wo = 7 for short relaxation times. This is clearly due 
to mass segregation from dynamical friction. This limits the 
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accuracy of the approximation. The maximum absolute dif- 
ference in the retention fractions are 0.12 for x = 2 and 0.08 
for x = 1. 

Bearing in mind these systematic differences, we can ap- 
proximate a retention fraction curve for models unlike those 
given in Table ^ by using the appropriately scaled retention 
curves from Table The effect of this approximation on the 
final retention rates will depend on the kick velocity distri- 
bution employed. A feeling for the size of the uncertainty 
can be taken from comparisons of rates based on Table U 
curves with rates based on scaled Table |l] curves. 



3 USE OF THE RETENTION FRACTIONS 

Once /fet(ffe) is available, we can calculate the retention 
fraction of neutron stars in a given model by assuming a 
distribution of kick velocities. Here, I use the distribution 
estimated by Lyne & Lorimer (1994) for the Population I 



pulsars. Note that the correct value of the parameter m is 
0.3, not 0.13 as published (D. Lorimer, private communi- 
cation). The results for a range of models are listed in the 
column labelled f" in Table q. For the most part, less than 
1% of the neutron stars are retained by the models. Only 
models with initial masses of 5 x 1O 6 M0 retain more than 
3% of their initial neutron stars. In no case does a model 
retain more than 10%. I have also tried models with a larger 
range of neutron star progenitors. This only changes the re- 
tention fractions by 10% or so, but such a model would have 
a larger number of neutron stars to retain. 

The number of neutron stars retained for each model, 
under the assumption that all the progenitors were single 
stars, is listed in the column labelled N^s m Table |(| Only 
models with x = 1, M — 5 x 1O 6 M and either highly 
concentrated (high Wo) or compact (small n) retain more 
than 1000 neutron stars in this scenario. For less massive 
models, the maximum number of retained neutron stars is 
of order 100 for x = 1 and of order 10 for x — 2. Clearly, 
under these assumptions, very few neutron stars would be 
retained in most globular clusters. 

Until now, I have assumed that the neutron star pro- 
genitors are single stars. Most massive stars, however, are 
members of multiple star systems, and the effects of multi- 
plicity on the neutron star retention rates must be included. 
There are three possible outcomes when the primary star in 
the binary becomes a supernova and receives its kick. The 
system may become unbound, the system may remain bound 
with a change in its center-of-mass velocity, or the neutron 
star may pass close to the companion, or even hit it in the 
extreme case, triggering mass transfer and probable merger, 
f ormin g a Thorne-Zytkow object. Brand t and Podsiadlowski 
(1995) used the same Lyne & Lorimer (1994) velocity spec- 



estimates given here are intended as suggestive rather than 
definitive. 

Alt hough not explicitly given in Brandt and Podsiad- 
lowski (1995), the same calculation can also give the re- 
sulting velocities for each of the components of a disrupted 
system and for the merged remnant of the TZO. The veloc- 
ity distributions f or the four types of results, assuming the 
Lyne & Lorimer (1994) kick distribution, are show in Fig- 
ure ^| The upper panel shows the distribution of speeds for 
the low velocity products: the former companion in the un- 
bound systems, the centre-of-mass motion of the bound bi- 
naries, and the motion of the merged objects. The bi-modal 
distribution of speeds for the unbound companion is a result 
of the two possible masses of the star under the two differ- 
ent regimes of stable or dynamical mass transfer (see Brandt 
and Podsiadlowski 1995 for details). The velocity distribu- 
tion for the new-born, unbound neutron stars is given in the 
l ower panel. It has the same peak as the Lyne & Lorimer 
( 1994 ) velocities, but a longer high- velocity tail. In passing, 
this may suggest that the high-velocity pulsars in the Lyne 
& Lorimer sample may have their origin in binaries and that 
the true kick-velocity distribution truncates at lower veloci- 
ties. If this were the case, the neutron star retention fraction 
would be somewhat increased, but since the mode of the dis- 
tributions remains at around 250 km s~ 2 , it would not be a 
large effect. 

For each channel of evolution, we can calculate the over- 
all probability of retention using the curves as from Table ^. 
For a disrupted system, the overall retention fraction is cal- 
culated using the velocity distributions for the freed neutron 
star and its former companion. The companion, if retained 
by the cluster, is assumed, by two-body interactions, to re- 
join the general background distribution of velocities. (Dy- 
namical friction is a quick process.) Its retention probability 
in its subsequent supernova is given by the single star reten- 
tion rate. The retention fraction of stars in the TZO channel 
is also calculated using their distribution of resultant veloc- 
ities. 

For the systems which remain bound, the distribution 
of changes in c enter -of-mass velocity based on Brandt and 
Podsiadlowski ( 1995 ) is used for the distribution of kick ve- 
locities. The subsequent evolution of these systems is highly 
dependent on their new orbital parameters and the environ- 
ment, but somewhere between and 2 neutron stars remain 
bound to the cluster from each system. Systems with suf- 
ficiently short periods are expected to become high-mass 
X-ray binaries and experience a spiral-in phase resulting in 
the formation of a TZO. The envelope will then be ejected 
leaving a single neutron star bound to the cluster, no sec- 
ond supernova having occured. Systems with longer periods 
will probably suffer a common-envelope phase and leave a 
short-period binary consisting of a helium star and a neutron 



trum to calculate the relative proportions of each outcome star. Such systems would have a second supernova (Bhat 



and the velocity distribution of the resulting binaries. 



tacharya & van den Heuvel 1991). The critical period is 



I h ave used the results for the high-mass x-ray binaries somewhere around 100 days 



Bodenheimer, & Os- 



from section 3.2 of Brandt and Podsiadlowski (1995). This triker 197S) and the majority of the systems in Brandt & 



work was meant as an exploratory investigation, so they used 
only one set of initial masses, but a distribution of initial 
periods. The calculations were based on a set of plausible, 
but arbitrary assumptions with regard to mass transfer in 
the evolving binary. As in Brandt and Podsiadlowski, the 



Podsiadlowski (1995) had periods less than 100 days. On the 



other hand, their Figure 8b shows that there is a correlation 
between centre-of-mass velocity and period in the sense that 
the longer period systems had lower centre-of-mass veloci- 
ties and thus were more likely to be bound to the cluster. 
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Overall, it is plausible that each initial binary results in a 
single neutron star remaining in the cluster. 

Given the fractions from each channel, the total fraction 
of neutron stars originating in binaries which are retained 
by the cluster is given by 

/ r bi » = (1 - ,f 6 )(/," 1 + frff) + Mftfr + (1 - ft)af b r ), (11) 

where ft, is the fraction of systems which remain bound, ft 
is the fraction of the bound systems becoming TZOs (the 
fraction of unbound systems which suffer close encounters 
as the neutron star escapes is very small), and the factors 
subscripted r are the retention fractions for the various pop- 
ulations: fr 1 for the unbound neutron stars, /" 2 for their for- 
mer companions, fr for the single neutron stars as above, f r 
for the TZOs, and fr for the systems which remain bound. 
The factor a deals with the fate of the binary systems. Its 
value lies between and 2, but is likely to be 1 as discussed 
above. For the models listed in Table I list in Table [] the 
retention fractions for the various processes. I have taken 
fb — .27 and ft = 0.26 based on Brandt & Podsiadlowski 
(1995), but it should be noted that these came from a partic- 
ular set of assumptions and initial conditions and should be 
considered as approximate. The first column gives the model 
number as in Table |^, the next five columns give fraction of 
neutron stars retained through each individual channel listed 
above. The next two columns give the combined fraction re- 
tained from the binary systems using eq.(pT|) and assuming 
that a = or 2. The next column gives the number of neu- 
tron stars retained if they are all single stars and the final 
three columns give the expected number of retained neutron 
stars under the assumption that all the neutron star progen- 
itors are in binary systems, but that either the bound bina- 
ries leave no neutron stars in the cluster (the a = case), 
iVrvs(O), that they leave one neutron star (the a = 1 case), 
N%f S (l), or that they leave two neutron stars in the cluster 
(the a = 2 case), N^, s (2). The a = 1 case gives the most 
likely value. 

By comparing the f" column and the two /, bln columns, 
it is clear that neutron stars originating from binaries are 
more likely to be retained than neutron stars from single 
stars. Further, the major mechanism for retaining neutron 
stars is through the channel of binaries which remain bound 
and binary after the first supernova. If the fate of these 
systems is as discussed above, then these binaries will ul- 
timately leave one neutron star without any further violent 
events to eject it from the cluster. Thus, on the assumptions 
made here, most globular cluster neutron stars will have 
passed through a high-mass X-ray binary phase. 

Between clusters, it is also clear from Table ^ that 
the shape of the initial mass function will play a crucial 
role in determining the number of neutron stars in the 
cluster. The variation in the retention fractions with the 
mass function slope is relatively small, so of two clusters 
with the same structural parameters — initial mass, size, and 
concentration — the cluster with the larger number of mas- 
sive stars will retain the larger number of neutron stars. This 
holds even if this cluster is somewhat less massive, less con- 
centrated, or larger than the other cluster with the steeper 
initial mass function. 

To illustrate this point further, and as an example of the 
use of the Michie-King model approximation, I estimate the 
number of retained neutron stars for two globular clusters. 



The main problem in making such estimates for particular 
globular clusters is that the initial conditions of the cluster 
are required, not the present state. Since the initial condi- 
tions are unknown, this presents a problem, but estimates 
can be made. I will consider two globular clusters: NGC 6397 
and u) Cen. For NGC 6397, approximate initia l conditions 
are known from the modelling in Drukier (1995). The model 



used is model U30-B. Since the cluster had a fairly low ini- 
tial mass, it is to be expect ed th at it lost a large majority of 
its neutron stars. Drukier ( 1995 ) showed that some neutron 
stars had to have been retained, so I have taken the model 
with the largest initial number of massive stars. The u) Cen 
model assumes that the cluster hasn't evolved very much, 
which is reasonable enough given its high mass and low con- 
centration, but does include the appropriate change in the 
total mass. The parameters I adopt for the initial models 
are listed in Table [?]. Table [?] also contains the estimates of 
the original and retained numbers of neutron stars under the 
same assumptions as in Table |(| except that for the number 
originating from binaries, a = 1 is assumed. NGC 6397's 
lower initial mass is balanced by its higher concentration, 
smaller size and larger number of neutron-star progenitors. 
ui Cen retains a much higher fraction of its neutron stars, 
but the two clusters end up retaining under the various sce- 
narios the same numbers of neutron stars. 



4 CONCLUSIONS 

I have presented tables based on Fokker-Planck models 
which give the retention fractions of neutron stars as a func- 
tion of the initial parameters of the models and of the neu- 
tron star kick velocity. I have also introduced an approxi- 
mation technique using Michie-King models. These can be 
combined with a distribution of kick velocities to estimate 
the number of neutron stars retained by a globular cluster. 
One difficulty with this approach, is that the present struc- 
ture of a globular cluster is not necessarily a good estimate of 
the state of the cluster when the neutron stars were formed. 

These retention fractions have been integra ted ov er the 
kick velocity distribution of Lyne & Lorimer (1994), and, 
as expected, very few single neutron stars are retained. The 
same distribution has been applied to the evolution of a 
model sample of high-mass binaries and the fraction of neu- 
tron stars retained from binaries was similarly estimated. 
Based on these estimates, the dominant source of globular 
cluster neutron stars would be via high-mass X-ray binaries. 
Both for neutron stars originating as single stars and in bi- 
naries, the number retained is strongly dependent on the 
number of neutron-star progenitors. Increasing the number 
of progenitors increases the total mass loss and reduces the 
fraction retained, but the larger number of stars more than 
makes up for this reductions. 

The retention fractions given in this paper contribute 
only to the first stage in estimating the observed population 
of millisecond pulsars and LMXBs in globular clusters. The 
details of the rest of the evolution; i.e. how the neutron stars 
are spun up, is a matter for further study, but this paper 
should provide the tools to estimate the number of neutron 
stars the clusters have to work with. 
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Table 2. Initial parameters for Fokkcr-Planck models 
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Table 4. Retention fractions for Fokker-Planck models in Table 
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1.000 


1.000 


1.000 


1 


1.000 


1.000 


1.000 


1. 


.000 


1.000 


1.000 


1.000 


2 


1.000 


1.000 


0.998 


1. 


.000 


1.000 


1.000 


0.998 


3 


0.998 


1.000 


0.984 


0. 


.995 


0.998 


0.995 


0.986 


4 


0.989 


0.997 


0.949 





.979 


0.988 


0.980 


0.953 


5 


0.964 


0.985 


0.889 


0. 


.944 


0.964 


0.945 


0.896 


6 


0.920 


0.964 


0.810 


0. 


.888 


0.932 


0.891 


0.820 


7 


0.850 


0.930 


0.721 





.809 


0.879 


0.812 


0.732 


8 


0.758 


0.874 


0.626 


0. 


.714 


0.812 


0.718 


0.640 


9 


0.650 


0.805 


0.534 


0. 


.609 


0.736 


0.615 


0.548 


10 


0.533 


0.719 


0.443 


0. 


.501 


0.648 


0.507 


0.458 


11 


0.417 


0.625 


0.356 


0. 


.394 


0.557 


0.401 


0.372 


12 


0.309 


0.521 


0.275 





.295 


0.468 


0.302 


0.292 


13 


0.217 


0.416 


0.201 


0. 


.207 


0.376 


0.214 


0.217 


11 


0.146 


0.312 


0.138 


0. 


.139 


0.283 


0.144 


0.153 


15 


0.094 


0.220 


0.089 





.089 


0.193 


0.093 


0.103 


16 


0.058 


0.143 


0.055 


0. 


.054 


0.117 


0.057 


0.065 


17 


0.035 


0.089 


0.032 


0. 


.032 


0.067 


0.034 


0.040 


18 


0.020 


0.053 


0.018 


0. 


.018 


0.037 


0.020 


0.023 


19 


0.011 


0.031 


0.010 





.010 


0.021 


0.011 


0.013 


20 


0.006 


0.017 


0.005 





.005 


0.011 


0.006 


0.007 


21 


0.003 


0.009 


0.003 


0. 


.003 


0.005 


0.003 


0.004 


22 


0.001 


0.005 


0.001 


0. 


.001 


0.003 


0.001 


0.002 


23 


0.001 


0.002 


0.001 





.001 


0.001 


0.001 


0.001 


24 


0.000 


0.001 


0.000 


0. 


.000 


0.000 


0.000 


0.000 


25 


0.000 


0.000 


0.000 


0. 


.000 


0.000 


0.000 


0.000 


26 


0.000 


0.000 


0.000 


0. 


.000 


0.000 


0.000 


0.000 



a Multiply index by this to get the in 





8 




9 


10 


11 


12 


13 


14 


15 


3.84 


3.84 


1.28 


2.56 


2.24 


5.52 


3.68 


1.2 


1 


.000 


1 


.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1 


.000 


1 


.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 





.998 


1 


.000 


1.000 


0.998 


1.000 


0.997 


0.998 


0.998 





.986 





.996 


0.998 


0.984 


0.999 


0.983 


0.983 


0.987 





.959 





.981 


0.989 


0.949 


0.990 


0.947 


0.945 


0.956 





.905 





.947 


0.973 


0.889 


0.966 


0.894 


0.881 


0.907 





.834 





.893 


0.940 


0.809 


0.923 


0.817 


0.798 


0.836 





.756 





.818 


0.894 


0.720 


0.858 


0.733 


0.708 


0.756 





.672 





.725 


0.832 


0.628 


0.768 


0.647 


0.613 


0.676 





.590 





.624 


0.762 


0.536 


0.662 


0.562 


0.519 


0.593 





.514 





.516 


0.683 


0.447 


0.547 


0.482 


0.428 


0.514 





.439 





.410 


0.605 


0.360 


0.431 


0.407 


0.340 


0.439 





.374 





.309 


0.520 


0.277 


0.323 


0.342 


0.258 


0.374 





.314 





.218 


0.435 


0.200 


0.228 


0.282 


0.183 


0.314 





.261 





.146 


0.340 


0.136 


0.155 


0.228 


0.122 


0.258 





.206 





.093 


0.237 


0.086 


0.100 


0.171 


0.077 


0.203 





.134 





.056 


0.137 


0.052 


0.063 


0.101 


0.046 


0.139 





.065 





.033 


0.073 


0.030 


0.038 


0.049 


0.026 


0.072 





.030 





.019 


0.038 


0.016 


0.022 


0.024 


0.014 


0.035 





.016 





.010 


0.019 


0.008 


0.012 


0.013 


0.007 


0.018 





.009 





.005 


0.010 


0.004 


0.006 


0.007 


0.004 


0.010 





.005 





.003 


0.005 


0.002 


0.003 


0.004 


0.002 


0.005 





.002 





.001 


0.002 


0.001 


0.002 


0.002 


0.001 


0.003 





.001 





.001 


0.001 


0.000 


0.001 


0.001 


0.000 


0.002 





.000 





.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.001 





.000 





.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 




23 




24 


25 


26 


27 


28 


29 


30 


12.32 




4 


g 


2.72 


2.76 


1.6 


7.2 


2. 1 


1 


.000 


1 


.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1 


.000 


1 


.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 





.998 


1 


.000 


0.998 


0.999 


0.999 


0.999 


1.000 


1.000 





.985 





.998 


0.984 


0.989 


0.995 


0.991 


0.998 


0.999 





.951 





.990 


0.949 


0.958 


0.979 


0.968 


0.989 


0.995 





.897 





.968 


0.891 


0.902 


0.941 


0.925 


0.963 


0.979 





.823 





.934 


0.811 


0.834 


0.882 


0.864 


0.917 


0.956 





.741 





.882 


0.723 


0.753 


0.801 


0.788 


0.843 


0.903 





.656 





.813 


0.629 


0.671 


0.704 


0.708 


0.749 


0.840 





.574 





.740 


0.536 


0.585 


0.597 


0.625 


0.639 


0.759 





.494 





.656 


0.445 


0.507 


0.488 


0.541 


0.521 


0.660 





.420 





.566 


0.359 


0.433 


0.380 


0.459 


0.405 


0.555 





.352 





.474 


0.277 


0.367 


0.282 


0.381 


0.299 


0.445 





.292 





.382 


0.204 


0.306 


0.197 


0.307 


0.209 


0.345 





.234 





.292 


0.142 


0.245 


0.131 


0.237 


0.140 


0.247 





.170 





.205 


0.093 


0.186 


0.083 


0.177 


0.090 


0.168 





.103 





.132 


0.058 


0.118 


0.050 


0.126 


0.056 


0.106 





.052 





.079 


0.035 


0.061 


0.029 


0.086 


0.033 


0.067 





.027 





.045 


0.020 


0.032 


0.017 


0.056 


0.019 


0.040 





.013 





.026 


0.011 


0.015 


0.009 


0.035 


0.010 


0.023 





.007 





.015 


0.006 


0.008 


0.005 


0.022 


0.005 


0.012 





.004 





.007 


0.003 


0.004 


0.002 


0.013 


0.003 


0.006 





.002 





.004 


0.002 


0.002 


0.001 


0.008 


0.001 


0.003 





.001 





.002 


0.001 


0.001 


0.001 


0.004 


0.001 


0.001 





.000 





.001 


0.000 


0.001 


0.000 


0.002 


0.000 


0.001 





.000 





.000 


0.000 


0.000 


0.000 


0.001 


0.000 


0.000 





.000 





.000 


0.000 


0.000 


0.000 


0.001 


0.000 


0.000 


km 


s" 1 


corresponding to the table entry. 
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Table 4. - continued. 



Model 


31 


32 


33 


34 


35 


36 


37 


v a 


1.6 


3.84 


1.56 


2.48 


3.36 


2.16 


6.72 


index 
















o 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


2 


1.000 


0.999 


1.000 


0.998 


1.000 


1.000 


0.998 


3 


0.998 


0.988 


0.995 


0.983 


1.000 


0.999 


0.983 


4 


0.989 


0.965 


0.980 


0.946 


0.994 


0.991 


0.945 


5 


0.964 


0.913 


0.946 


0.885 


0.979 


0.971 


0.880 





0.917 


0.851 


0.890 


0.802 


0.951 


0.932 


0.795 


7 


0.847 


0.774 


0.813 


0.710 


0.909 


0.871 


0.702 


8 


0.755 


0.691 


0.719 


0.614 


0.844 


0.789 


0.603 


9 


0.645 


0.610 


0.616 


0.519 


0.763 


0.690 


0.507 


10 


0.527 


0.533 


0.509 


0.426 


0.666 


0.580 


0.415 


11 


0.412 


0.457 


0.402 


0.340 


0.559 


0.468 


0.327 


12 


0.306 


0.388 


0.304 


0.258 


0.453 


0.361 


0.247 


13 


0.216 


0.325 


0.217 


0.187 


0.349 


0.266 


0.176 


14 


0.145 


0.257 


0.148 


0.128 


0.254 


0.187 


0.119 


15 


0.094 


0.193 


0.096 


0.083 


0.174 


0.127 


0.076 


16 


0.059 


0.125 


0.060 


0.051 


0.114 


0.083 


0.046 


17 


0.035 


0.066 


0.036 


0.030 


0.072 


0.053 


0.027 


18 


0.020 


0.033 


0.021 


0.017 


0.044 


0.032 


0.016 


19 


0.011 


0.017 


0.012 


0.010 


0.025 


0.019 


0.008 


20 


0.006 


0.009 


0.006 


0.005 


0.014 


0.011 


0.005 


21 


0.003 


0.005 


0.003 


0.003 


0.008 


0.006 


0.002 


22 


0.001 


0.002 


0.002 


0.001 


0.004 


0.003 


0.001 


23 


0.001 


0.001 


0.001 


0.001 


0.002 


0.002 


0.001 


24 


0.000 


0.001 


0.000 


0.000 


0.001 


0.001 


0.000 


25 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 



Model 


46 


47 


48 


49 


50 


51 


52 


v a 


8.64 


2.88 


5.76 


2 


1.84 


1.68 


7.44 


index 



















1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


2 


0.999 


1.000 


0.998 


0.998 


1.000 


1.000 


1.000 


3 


0.990 


0.998 


0.984 


0.984 


0.998 


0.997 


0.999 


4 


0.969 


0.986 


0.945 


0.949 


0.989 


0.987 


0.995 


5 


0.927 


0.963 


0.884 


0.887 


0.963 


0.963 


0.980 


6 


0.860 


0.931 


0.800 


0.807 


0.919 


0.922 


0.954 


7 


0.788 


0.868 


0.708 


0.718 


0.846 


0.861 


0.912 


8 


0.709 


0.790 


0.611 


0.622 


0.755 


0.785 


0.842 


9 


0.628 


0.708 


0.515 


0.527 


0.647 


0.696 


0.761 


10 


0.544 


0.617 


0.423 


0.434 


0.531 


0.600 


0.665 


11 


0.463 


0.524 


0.335 


0.347 


0.415 


0.502 


0.557 


12 


0.387 


0.429 


0.254 


0.267 


0.310 


0.405 


0.451 


13 


0.314 


0.335 


0.183 


0.195 


0.219 


0.314 


0.350 


14 


0.245 


0.246 


0.126 


0.136 


0.148 


0.234 


0.259 


15 


0.179 


0.168 


0.082 


0.089 


0.097 


0.167 


0.180 


16 


0.113 


0.107 


0.051 


0.056 


0.061 


0.114 


0.117 


17 


0.061 


0.064 


0.030 


0.034 


0.037 


0.076 


0.076 


18 


0.033 


0.037 


0.018 


0.020 


0.022 


0.049 


0.047 


19 


0.017 


0.022 


0.010 


0.011 


0.012 


0.031 


0.028 


20 


0.009 


0.011 


0.005 


0.006 


0.006 


0.018 


0.016 


21 


0.005 


0.006 


0.003 


0.003 


0.003 


0.011 


0.009 


22 


0.002 


0.003 


0.001 


0.002 


0.002 


0.006 


0.004 


23 


0.001 


0.002 


0.001 


0.001 


0.001 


0.003 


0.002 


24 


0.001 


0.001 


0.000 


0.000 


0.000 


0.002 


0.001 


25 


0.000 


0.000 


0.000 


0.000 


0.000 


0.001 


0.000 


26 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 



a Multiply index by this to get the in 



38 


39 


40 


41 




42 




43 




44 


45 


2.16 


2.16 


9 


6.16 


2.04 


1.36 


1.28 


3.04 


1.000 


1.000 


1.000 


1.000 


1 


000 


1 


.000 


1 


.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1 


000 


1 


.000 


1 


.000 


1.000 


0.999 


1.000 


1.000 


1.000 


1 


000 


1 


.000 


1 


.000 


0.999 


0.991 


0.996 


0.998 


0.995 





.998 





.995 





.999 


0.993 


0.968 


0.982 


0.988 


0.977 





.988 





.980 





.990 


0.968 


0.931 


0.951 


0.969 


0.942 





.963 





.944 





.967 


0.932 


0.870 


0.898 


0.930 


0.882 





.927 





.887 





.924 


0.879 


0.798 


0.824 


0.877 


0.801 


0. 


869 





.808 





.858 


0.810 


0.719 


0.736 


0.804 


0.702 





.793 





.713 





.769 


0.732 


0.641 


0.635 


0.723 


0.596 





712 





.608 





.662 


0.655 


0.560 


0.530 


0.631 


0.486 





.617 





.499 





.550 


0.572 


0.487 


0.425 


0.538 


0.379 





.524 





.393 





.436 


0.497 


0.414 


0.326 


0.444 


0.280 





427 





.295 





.330 


0.422 


0.346 


0.239 


0.353 


0.197 





334 





.210 





.236 


0.347 


0.280 


0.166 


0.263 


0.131 





.246 





142 





.162 


0.282 


0.215 


0.110 


0.182 


0.084 





.168 





.092 





108 


0.214 


0.150 


0.071 


0.117 


0.051 





.106 





.057 





.069 


0.152 


0.088 


0.044 


0.070 


0.030 





064 





.034 





.042 


0.092 


0.048 


0.026 


0.043 


0.017 





038 





.020 





.025 


0.050 


0.025 


0.015 


0.024 


0.010 





.022 





Oil 





.014 


0.028 


0.014 


0.008 


0.014 


0.005 





Oil 





.006 





.008 


0.015 


0.007 


0.005 


0.007 


0.002 





006 





.003 





.004 


0.008 


0.004 


0.002 


0.004 


0.001 





003 





.001 





.002 


0.004 


0.002 


0.001 


0.002 


0.001 





001 





.001 





.001 


0.002 


0.001 


0.000 


0.001 


0.000 





.001 





.000 





.000 


0.001 


0.000 


0.000 


0.000 


0.000 





000 





.000 





.000 


0.000 


53 


54 


55 


56 




57 




58 




59 


60 


5 


1.6 


4.96 


1.6 


2.64 


2.4 


6.96 


2.28 


1.000 


1.000 


1.000 


1.000 


1 


000 


1 


.000 


1 


.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1 


000 


1 


.000 


1 


.000 


1.000 


1.000 


1.000 


1.000 


1.000 





.999 


1 


.000 





.999 


1.000 


0.998 


1.000 


0.995 


0.998 





.992 


1 


.000 





.990 


0.998 


0.990 


0.994 


0.980 


0.990 





.963 





.994 





.970 


0.987 


0.965 


0.983 


0.945 


0.968 





.921 





.978 





.923 


0.965 


0.923 


0.957 


0.887 


0.935 





858 





.944 





.861 


0.932 


0.851 


0.920 


0.809 


0.882 





778 





.887 





.788 


0.876 


0.761 


0.861 


0.712 


0.813 





.697 





.821 





.709 


0.804 


0.654 


0.783 


0.607 


0.729 





616 





.734 





.626 


0.719 


0.539 


0.695 


0.500 


0.639 





.529 





.631 





.540 


0.632 


0.425 


0.597 


0.392 


0.551 





451 





.526 





.458 


0.535 


0.319 


0.496 


0.295 


0.457 





.374 





.417 





.378 


0.442 


0.228 


0.398 


0.210 


0.368 





301 





.317 





.307 


0.350 


0.155 


0.304 


0.142 


0.280 





.231 





.228 





.236 


0.266 


0.102 


0.221 


0.092 


0.202 





165 





155 





.173 


0.190 


0.065 


0.151 


0.058 


0.136 


0. 


105 





101 





.114 


0.127 


0.040 


0.101 


0.035 


0.086 


0. 


059 





.065 





.069 


0.082 


0.023 


0.065 


0.020 


0.054 





.033 





.039 





.039 


0.050 


0.013 


0.041 


0.011 


0.033 





017 





.023 





.021 


0.030 


0.007 


0.025 


0.006 


0.018 





009 





.013 





.011 


0.018 


0.004 


0.014 


0.003 


0.010 





.005 





.006 





.006 


0.009 


0.002 


0.008 


0.001 


0.005 





002 





.003 





.003 


0.005 


0.001 


0.004 


0.001 


0.003 





001 





.001 





.002 


0.003 


0.000 


0.002 


0.000 


0.001 





000 





.000 





.001 


0.001 


0.000 


0.001 


0.000 


0.001 





000 





.000 





.000 


0.001 


0.000 


0.000 


0.000 


0.000 





.000 





.000 





.000 


0.000 



km s 1 corresponding to the table entry. 
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Table 5. Central escape velocity as a function of Wo 



Wn 

VV (J 


Vfi (km s ^ 




1 


0.18 


0.65 


2 


0.19 


0.72 


Q 
O 


n 9n 
u.zu 


U.ol 


4 


0.22 


0.93 


5 


0.25 


1.09 


6 


0.29 


1.29 


7 


0.35 


1.55 


8 


0.41 


1.85 


9 


0.46 


2.12 


10 


0.47 


2.36 


11 


0.47 


2.55 
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Table 6. Retention fractions and retained numbers of neutron stars for the models in Table 



Run 


fr 


fr 


ft 


fr 1 


fr 2 


/r in (0) 


/r in ( 2 ) 


N s „ 

N S 


IV b K ' 


IV o v ' 


IV o K > 


1 


0.002 


0.041 


0.002 


0.001 


0.048 


0.001 


0.009 


3 


2 


9 


15 


2 


0.009 


0.306 


0.057 


0.006 


0.357 


0.010 


0.072 


77 


86 


350 


620 


3 


0.001 


0.014 


0.000 


0.000 


0.010 


0.000 


0.003 


2 





3 


5 


4 


0.017 


0.532 


0.187 


0.013 


0.555 


0.029 


0.135 


290 


500 


1400 


2300 


5 


0.003 


0.066 


0.005 


0.001 


0.089 


0.001 


0.015 








1 


1 


6 


0.006 


0.183 


0.023 


0.003 


0.244 


0.005 


0.041 


51 


43 


200 


350 


7 


0.078 


0.914 


0.754 


0.061 


0.905 


0.149 


0.332 


6700 


13000 


21000 


28000 


8 


0.012 


0.395 


0.099 


0.008 


0.433 


0.016 


0.095 


5 


6 


23 


38 


9 


0.011 


0.372 


0.078 


0.007 


0.425 


0.014 


0.088 


190 


240 


870 


1500 


10 


0.002 


0.032 


0.001 


0.001 


0.031 


0.001 


0.007 











1 


11 


0.005 


0.141 


0.017 


0.002 


0.189 


0.004 


0.032 


13 


34 


150 


270 


12 


0.004 


0.121 


0.012 


0.002 


0.167 


0.003 


0.027 


34 


26 


130 


230 


13 


0.022 


0.606 


0.282 


0.017 


0.620 


0.042 


0.163 


18 


34 


82 


130 


11 


0.009 


0.298 


0.056 


0.005 


0.349 


0.010 


0.070 


150 


170 


680 


1200 


15 


0.001 


0.019 


0.001 


0.000 


0.018 


0.000 


0.004 














16 


0.008 


0.267 


0.042 


0.005 


0.335 


0.008 


0.062 


140 


140 


600 


1100 


17 


0.001 


0.013 


0.000 


0.000 


0.008 


0.000 


0.003 














18 


0.003 


0.084 


0.007 


0.001 


0.113 


0.002 


0.019 


26 


17 


86 


160 


19 


0.049 


0.876 


0.633 


0.040 


0.862 


0.104 


0.279 


4200 


8900 


16000 


24000 


20 


0.008 


0.267 


0.041 


0.005 


0.333 


0.008 


0.062 


3 


3 


14 


25 


21 


0.003 


0.075 


0.006 


0.001 


0.099 


0.002 


0.017 


26 


17 


77 


150 


22 


0.006 


0.194 


0.027 


0.003 


0.247 


0.005 


0.044 


100 


86 


430 


750 


2.3 


0.102 


0.936 


0.807 


0.079 


0.928 


0.183 


0.370 


410 


740 


1100 


1500 


24 


0.014 


0.496 


0.136 


0.010 


0.529 


0.023 


0.122 


11 


18 


58 


98 


25 


0.040 


0.789 


0.514 


0.032 


0.785 


0.082 


0.240 


3400 


7000 


14000 


21000 


26 


0.006 


0.199 


0.028 


0.003 


0.250 


0.006 


0.045 


2 


2 


10 


18 


27 


0.006 


0.180 


0.023 


0.003 


0.239 


0.005 


0.041 


100 


86 


390 


700 


28 


0.002 


0.055 


0.004 


0.001 


0.071 


0.001 


0.012 


17 


9 


60 


100 


29 


0.035 


0.815 


0.490 


0.029 


0.802 


0.076 


0.239 


3000 


6500 


13000 


20000 


30 


0.005 


0.177 


0.021 


0.003 


0.242 


0.005 


0.040 


2 


2 


9 


16 


31 


0.002 


0.044 


0.002 


0.001 


0.050 


0.001 


0.010 


17 


9 


43 


86 


32 


0.012 


0.404 


0.101 


0.008 


0.441 


0.017 


0.098 


10 


11 


16 


79 


33 


0.002 


0.039 


0.002 


0.001 


0.044 


0.001 


0.009 


17 


9 


43 


77 


34 


0.004 


0.127 


0.014 


0.002 


0.170 


0.003 


0.028 


68 


51 


270 


480 


35 


0.010 


0.364 


0.071 


0.007 


0.423 


0.013 


0.086 


8 


10 


39 


69 


36 


0.004 


0.121 


0.012 


0.002 


0.165 


0.003 


0.027 


68 


51 


260 


460 


37 


0.027 


0.679 


0.352 


0.021 


0.685 


0.053 


0.189 


2300 


4500 


10000 


16000 


38 


0.004 


0.130 


0.014 


0.002 


0.177 


0.003 


0.029 


2 


1 


6 


12 


39 


0.004 


0.110 


0.011 


0.002 


0.150 


0.002 


0.024 


68 


34 


220 


410 


10 


0.067 


0.928 


0.759 


0.053 


0.914 


0.137 


0.322 


270 


550 


920 


1300 


41 


0.025 


0.692 


0.326 


0.020 


0.694 


0.050 


0.188 


2100 


4300 


10000 


16000 


42 


0.004 


0.115 


0.011 


0.002 


0.159 


0.003 


0.025 


2 


1 


6 


10 


13 


0.002 


0.023 


0.001 


0.000 


0.021 


0.000 


0.005 


17 





26 


13 


11 


0.002 


0.020 


0.001 


0.000 


0.017 


0.000 


0.004 


17 





17 


34 


45 


0.008 


0.284 


0.049 


0.005 


0.337 


0.009 


0.066 


6 


7 


31 


53 


46 


0.056 


0.866 


0.653 


0.045 


0.858 


0.113 


0.286 


230 


450 


800 


1100 


17 


0.007 


0.254 


0.039 


0.004 


0.319 


0.008 


0.058 


6 


6 


27 


47 


48 


0.020 


0.588 


0.247 


0.015 


0.605 


0.038 


0.155 


1700 


3200 


8200 


13000 


49 


0.003 


0.075 


0.006 


0.001 


0.099 


0.002 


0.016 


51 


34 


150 


270 


50 


0.003 


0.069 


0.005 


0.001 


0.089 


0.001 


0.015 


51 


17 


140 


260 


51 


0.003 


0.067 


0.005 


0.001 


0.086 


0.001 


0.015 


51 


17 


110 


260 


52 


0.047 


0.888 


0.637 


0.038 


0.871 


0.103 


0.280 


190 


410 


770 


1100 


53 


0.018 


0.582 


0.201 


0.014 


0.603 


0.032 


0.148 


1500 


2700 


7700 


13000 


54 


0.003 


0.070 


0.005 


0.001 


0.090 


0.001 


0.015 


1 





3 


6 


55 


0.017 


0.547 


0.182 


0.013 


0.573 


0.029 


0.138 


1500 


2500 


7200 


12000 


56 


0.003 


0.064 


0.004 


0.001 


0.082 


0.001 


0.014 


1 





3 


6 


57 


0.006 


0.189 


0.025 


0.003 


0.243 


0.005 


0.043 


5 


■1 


19 


35 


58 


0.005 


0.170 


0.020 


0.003 


0.232 


0.004 


0.038 


4 


3 


17 


31 


59 


0.037 


0.775 


0.491 


0.029 


0.772 


0.077 


0.231 


150 


310 


620 


930 


60 


0.005 


0.157 


0.018 


0.003 


0.213 


0.004 


0.035 


4 


3 


16 


28 
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Table 7. Parameters and results for models of two clusters 

Cluster Wo r,(pc) x M (M ) / M iV Massive # / r bin (0) / r bin (2) iV^ s N b NS (l) 

NGC 6397 6 22 0.9 5.5 x 10 5 0.77 9700 0.005 0.005 0.074 52 380 

u Cen 5 67 2. 4.3 x 10 6 0.99 3500 0.012 0.018 0.193 42 370 
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Figure 1. Fraction of neutron stars retained as a function of relative kick velocity for a range of Michic-King models. From right to left 
the dimensionless central potentials are from Wo = 1 to Wo = 11 in unit increments. 
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Figure 2. These figures show the effect on the retention fractions of varying the initial model parameters. All the panels show the 
retention fraction as a function of kick velocity for a range of models differing in only one parameter. Solid lines are for x = 1 and dotted 
lines are for x = 2. Note that at a given v k , the x = 2 models retain a higher fraction than do the x = 1 models, (a) The models all 
have M = 10 6 Mq and Wo = 5. The retention fraction decreases with increasing model size as defined by r;. In order of increasing 
at constant / ro t, r ; = 97, 65, 45, 30, and 15 pc. (b) The models all have r ; = 30 pc and Wo = 5. The retention fraction increases with 
increasing mass. In order of increasing Vk at constant / re t, M = 10°, 5 X 10 , 10 6 , and 5 X 10 Mq. (c) The models all have r\ = 15 pc 
and M = 10 6 Mq. The retention fraction increases with increasing concentration. In order of increasing at constant / re t, Wo = 3, 5, 
and 7. 
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Figure 3. Distribution of final speeds for the post-supernova products of binary systems. Top: (dashed) Ex-companions from unbound 
systems, (dotted) Centre-of-mass motion of bound binaries, (solid) merged stars (Thorne-Zytkow objects). Bottom: Neutron stars from 
unbound systems. The dotted line is the 3D velocity distribution from Lyne & Lorimer (1994). 
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